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Recent advances in non-equilibrium statistical mechanics and single molecule 
technologies make it possible to extract free energy differences from irre- 
versible work measurements in pulling experiments. To date, free energy 
recovery has been focused on native or equilibrium molecular states, whereas 
free energy measurements of kinetic states (i.e. finite lifetime states that are 
generated dynamically and are metastable) have remained unexplored. Ki- 
netic states can play an important role in various domains of physics, such as 
nanotechnology or condensed matter physics. In biophysics, there are many 
examples where they determine the fate of molecular reactions: protein and 
peptide-nucleic acid binding, specific cation binding, antigen-antibody in- 
teractions, transient states in enzymatic reactions or the formation of tran- 
sient intermediates and non-native structures in molecular folders. Here we 
demonstrate that it is possible to obtain free energies of kinetic states by ap- 
plying extended fluctuation relations. This is shown by using optical tweezers 
to mechanically unfold and refold DNA structures exhibiting intermediate 
and misfolded kinetic states. 

Kinetic states are observed under non-equilibrium conditions and have higher free 
energies than native states. Yet, they can be crucial, as shown by the role that misfolded 
proteins play in numerous severe diseases [I]. The measurement of the free energy of 
formation of kinetic states is therefore a central question in biophysics. Recent theo- 
retical developments known as fluctuation relations [131 121 IH G3 E] have been applied 
to extract free energy differences of equilibrium states from irreversible work measure- 
ments. Applications include the measurement of the free energy of formation of RNA 
and DNA hairpins [7]; the determination of the stability of native domains in proteins 
[8] ; the measurement of mechanical torque in rotary motors [9] ; the conversion of infor- 
mation into work in systems under feedback control |1UJ ; or the recovery of free energy 
landscapes from unidirectional work measurements |11| H^] . 

The characterization of kinetic states under non-equilibrium conditions remains a 
challenging problem. Here we use a recently introduced extended fluctuation relation 
(EFR) to extract free energies of kinetic states and thermodynamic branches using irre- 
versible work measurements [HHIO]. m the EFR, a kinetic state is a partially equilibrated 
region of configurational space, meaning that during a finite timescale the system is con- 
fined and thermalized within that region [15] . This is mathematically described by a 
Boltzmann-Gibbs distribution restricted to configurations contained in that region (Fig. 
la). 

Let A, B denote any two kinetic states and A a control parameter. We consider a 
forward (F) non-equilibrium process, where the system starts in partial equilibrium in 
A at Ao, and its time-reversed (R), where the partial equilibrium condition is required 
over B at Ai. In the F process A varies from Ao to Ai during a time r according to a 
predetermined protocol A(i). For the R process the time-reversed protocol A(r — t) is 
used. The EFR reads fTUl: 
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where AGab = Gb(Ai) — Ga(Ao) is the free energy difference between kinetic states B 
at Ai and A at A ; Pp^ B {W) (Pr"~ B '(—W)) denotes the partial work distribution for 
the F (R) process over the fraction of paths 4>p~^ B (4>^~ B ) starting in A (B) at Ao (Ai) 
and ending in B (A) at Ai (Ao); &b is the Boltzmann constant and T the temperature 
of the environment. 

We applied equation (fTl) to extract free energy differences of kinetic states from 
mechanical unfolding/folding experiments performed on DNA hairpins, which are model 
systems easy to design and synthesize. Their free energies of formation can be predicted 
using the nearest-neighbor (NN) model with the unified-oligonucleotide (UO) set of 
parameters [3] E| or with recently derived energies from unzipping experiments [18] 
(Methods, Supplementary Section SI). Molecules exhibiting two types of kinetic states 
were investigated (Fig. lb): molecules II and 12 have intermediate kinetic states on- 
pathway to the native state, and molecules Ml and M2 have misfolded kinetic states 
off-pathway to the native state. To establish the validity of our approach, we first show 
results for molecules II and Ml where free energies measured from the EFR applied to 
non-equilibrium pulling experiments can be compared with free energies obtained from 
equilibrium hopping experiments. The method is then applied to molecules 12 and M2 
where irreversibility or low signal-to-noise ratio in hopping traces preclude equilibrium 
based free energy measurements. 

The experimental setup is shown in Fig. lc |18|l2] (Methods). We steer up and down 
the position of the optical trap to mechanically unfold and refold the DNA hairpin, and 
measure the force acting on the hairpin as a function of the relative trap-pipette distance, 
which is the control parameter A [20] . We measure the work as the area below the force- 
distance curve (hereafter referred as FDC, inset of Fig. 2a) along many trajectories. 
Throughout this paper unfolding (folding) corresponds to the F (R) process. 

First we apply equation (TH) to hairpin II characterized by three conformational states 
(Fig. lb, Supplementary Section SI): native (N), intermediate (I) and unfolded (U). Ex- 
perimental hopping traces measured under equilibrium conditions and non-equilibrium 
FDCs exhibit three force branches corresponding to the three states (Fig. 2a, Supple- 
mentary Video VII and Sections S2, S3). Fig. 2b shows the partial work distributions 
measured from a collection of FDCs by taking Ao=0, where ^4=N, and Ai=55.6 nm, 
where the three states are observed (B=N, I or U). These partial work distributions sat- 
isfy equation ([I]) (Supplementary Section S4). Hysteresis effects are stronger for B=l, 
U than for i?=N, as the timescale related to the pulling protocol is typically shorter 
than the timescale for crossing the kinetic barrier separating two states. The accep- 
tance ratio method |21[ [22] applied to extract the free energy differences between states 
gives AGnu — AGnn> an d AGni lies 2 k^T above (Fig. 2c, Methods). Fig. 2d shows 
the reconstruction of the three thermodynamic branches by fixing Ao=0 and varying Ai 
between 45 and 65 nm. The vertical dashed-dotted line at A c =55.6 nm indicates the 
coexistence point of N and U. The full equilibrium free energy of the system, defined 
as AG = -£; B Tlog ( e - AG N N /fc B T + e -AG N i/fc B T + e -AG NV /k B T^ hag algo been measurec i 

(Fig. 2d, black line). The right inset in Fig. 2d shows the free energy of each state 
measured relative to AG. For A < A c (A > A c ) N (U) is the most stable state, while I 



is never the absolute free energy minimum for any A. The left inset in Fig. 2d shows 
the contribution of 4>f I 4>r to the measured free energies throughout the A range. 
Dropping this term or misidentifying states along the FDC leads to wrong free energy 
predictions (Supplementary Sections S5, S6). 

By subtracting the elastic contributions due to stretching the handles and the released 
ssDNA (Methods, Supplementary Section S7) we extract the free energies of formation 
of the different structures with respect to the random coil state at zero force. We get 
AGjJju = 55 ± 3 U-qT and AG^j = 30 ± 3 k^T, in agreement with free energy predictions 
and results from equilibrium-based hopping experiments (Table 1). 

Next we study hairpin Ml, which can fold into two unrelated structures (Fig. lb, 
Supplementary Section SI): the native (N) and the misfolded (M). Equilibrium hopping 
experiments exhibit very fast kinetics and two clearly separated hopping regions: in one 
region N coexists with an intermediate state on-pathway; in the other region M and U 
coexist with another intermediate (Supplementary Section 2). For simplicity reasons, we 
chose not to characterize these intermediate states. In non-equilibrium experiments, two 
FDC patterns are identified corresponding to the two structures (Fig. 3a, Supplementary 
Video VM1 and Section S3). In contrast to the unfolding/folding cycles that start and 
end in N, those that start and end in M show almost no hysteresis (Fig. 3b), indicating 
low kinetic barriers between M and U. Due to kinetic competition of loop formation 
between M and N, M has a basin of attraction larger than N during folding (~80% 
of folding trajectories end in M), lower thermodynamic stability and larger molecular 
extension at low forces (Fig. 2a, inset). 

In Fig. 3c we apply the acceptance ratio method to recover the free energy differences 
between U at Ai=130 nm and M (lower set of measurements) or N (upper set) at 
Ao=0 nm. By subtracting the handles and ssDNA contributions we extract the free 
energy of formation of each structure at zero force, obtaining AGj^y = 47 ± 2 k^T and 
AGjJjjj = 62 ±3 k^T. The distribution of free energies for different molecules and pulling 
speeds is shown in Fig. 3d. The difference between the average of both distributions is in 
agreement with predictions based on the NN model and with results from equilibrium- 
based hopping experiments (Table 1). 

Finally, to illustrate the power of the method we show the case of molecules 12 and 
M2, where it is not possible to recover free energies using equilibrium based methods 
(Supplementary Section S3). 

Hairpin 12 has two intermediate states on pathway, hereafter referred as I' and I" |23j 
(Supplementary Section SI). In pulling experiments four force branches, corresponding 
to states N, I', I" and U are distinguished (Fig. 4a, Supplementary Video VI2 and 
Sections S3, S4). In order to measure forward and reversed partial work distributions 
for the four states we pull back and forth the molecule between Ao=0 (where the molecule 
is in equilibrium at N) and Ai=183 nm (where the molecule is partially equilibrated at 
states N, I', I", U; Fig. 4a,b). This protocol is subtly different from the standard pulling 
experiments we did for the rest of molecules (II, Ml, M2) where the molecule is never 
in an intermediate state at initial and final values of A. Due to the larger hysteresis 
exhibited by this molecule (Supplementary Section S7), the standard protocol does not 



generate reverse trajectories that sample all four states for any value of A. 

In table 1 we show the values of the free energy of formation of the different kinetic 
states obtained with the EFR. The size of the error bars is comparable to the discrepancy 
between the free energy predictions using the NN model with the UO set of parameters 
[31 H] and unzipping data [18] . To evaluate the free energy branches of the different states 
(Fig. 4c) we could repeat the experiment for different final values of A and measure the 
corresponding AGnb(A), B=N, I', I" or U. For simplicity we use an extended version of 
the Jarzynski equality (EJE) obtained by multiplying the EFR with the reversed work 
distribution and integrating over the work (Supplementary Section S8), 

AG AB = -k B Tlog ( p^J - k B Tlog(e «** \ ( 2 ) 

Equation pi) only requires data from the F process and we apply it to pulling experiments 
recorded by setting extreme values of A (light curves in Fig. 4a). Similarly to the 
Jarzynski estimator [13], the EJE is strongly biased [2]. To estimate the magnitude of 
the bias we took the difference between the free energy AGnb obtained using equation 
pi) to the one obtained with the acceptance ratio method in pulling experiments where 
kinetic states are partially equilibrated at Ai (dark curves in Fig. 4a). Therefore, from 
the free energy branches obtained using equation p| we subtracted this estimated bias 
for each state (we assumed it to be equal for all values of Ai, Supplementary Section S8). 
In contrast to II, kinetic intermediates found in 12 become the most stable states in a 
given range of A. For low values of A, stability is determined by N, and as A increases 
stability shifts to F, I" and finally to U (Fig. 4c). 

Hairpin M2 can fold into one native structure (N) and two misfolded structures 
(M', M") following alternative folding pathways (Fig. 5a, Supplementary Section SI). 
Whereas it is easy to identify trajectories that fold into N (red/blue FDCs in Fig. 5a 
left, ~50% of trajectories), distinguishing trajectories that misfold into M' or M" is not 
straightforward. Careful inspection reveals two different patterns of unfolding curves 
that start at a misfolded state: either the molecule unfolds quasi-reversibly without 
intermediates (purple FDC in Fig. 5a middle, ~30% of trajectories), or it folds back 
to N before it unfolds (cyan FDC in Fig. 5a right, ~20% of trajectories) [llj. We 
interpret the former as trajectories following the M'— >-U pathway and the latter following 
the M"— >N— >\J pathway (Supplementary Video VM2). This is supported by two facts 
(Supplementary Section S5). First, M' consists of four small hairpins that confer low 
mechanical stability to the structure that gently unfolds under tension. Second, M" has 
a large stem in common with N (Fig. 5a, top) which is surrounded by two small hairpins 
with low mechanical stability. Once these two hairpins unfold around 9-10 pN, force 
remains low enough for the molecule to fold back to N before unfolding. Combining 
equation (Fy), the partial work distributions (Fig. 5b) and handles and ssDNA elastic 
contributions leads to the free energy values AG^ = 94 ± 2, AG^jj = 60 ± 3 and 
AGy HU = 70±3 k B T, in agreement with theoretical predictions (Table 1). Fig. 5c shows 
the reconstruction of the four thermodynamic branches relative to the full equilibrium 
free energy of the system, AG = ~k B T log Y^a=n,M'M",u e~ AG - Av ^ kBT , by fixing Ai=230 



nm and varying Ao between and 150 nm. 

Summarizing, we have shown how the EFR can be used to extract free energies of 
non-equilibrium kinetic structures in DNA hairpins exhibiting intermediate and mis- 
folded states. The method accurately works in far from equilibrium situations and when 
equilibrium experiments are insufficient to characterize non-native states. There are two 
main differences between the EFR in equation (fTl) and the Crooks relation |12j : the 
partial work distributions and the prefactor <j)p /4'r~ ■> which introduces the addi- 
tional correction —k-QT\og{4>p B /4>^ B ) into the Crooks estimation of the free energy 
difference between kinetic states. The omission of such correction yields wrong relative 
thermodynamic stabilities for the free energy branches of the different kinetic states 
(Supplementary Section S7) [10]. Moreover, for the case of misfolded structures that 
apparently unfold/misfold reversibly (Fig. 3a right and Fig. 5a middle), AGmu is 
not just equal to the measured reversible work during unfolding since since the term 
kftT log (fift^^ must be added (here </> I p _ShU = 1 as F processes always end at U). Al- 
though this correction is small for states M (Ml) and M" (M2) (~ 0.2 k B T and 1.2 k B T 
respectively), it is important in situations where </>/j <C 1, even if very low hysteresis 
is obtained between the F and R processes. For example, the neglection of a 1% mis- 
folding probability would underestimate by 4.5 k B T the free energy of formation of the 
misfolded state. 

The main limitation of the method is the identification of kinetic states from the mea- 
sured signal. In this regard, a combination of fluorescence techniques, such as FRET, 
with force measurements, and the application of advanced statistical methods (e.g. hid- 
den Markov models or Bayesian inference) might be very useful. Our methodology 
should find many applications that range from molecular biophysics to condensed matter 
physics. Any situation where equilibrium experiments are unpractical should be treat- 
able with different versions of equation (1). To start with, the method can be employed 
for measuring free energies of kinetic structures that appear in many molecular reac- 
tions, such as RNA, proteins, and many kinetic states related to intermolecular binding, 
or transient non-equilibrium states that are essential in polymerization reactions (e.g. 
ATP or ADP bound states in motor proteins). 

Methods 

Molecular synthesis 

The designed DNA molecules linked to 29 bp dsDNA-handles were synthesized as de- 
scribed in [2]. For the specific attachments to the DNA molecular construction we used 
streptavidin-coated polystyrene microspheres (1.87 fim, Spherotech, Libertyville, IL) 
and protein G microspheres (3.0-3.4 /iin; G. Kisker Gbr, Products for Biotechnologie, 
Steinfurt, Germany) coated with anti-digoxigenin polyclonal antibodies (Roche Applied 
Science, Spain). Attachment to the anti-digoxigenin microspheres was achieved first by 
incubating the beads with the tether DNA. The second attachment was achieved in the 
fluidics chamber and was accomplished by bringing a trapped anti-digoxigenin and an 



immobilized streptavidin microsphere close to each other. 

Bennett acceptance ratio method 

This method is used to estimate the free energy difference AGab between two states 
from non-equilibrium work measurements. Given a set of npinii) forward(reversed) work 
measurements, it is shown in |21|, [22] that the solution of the following transcendental 
equation: 



f3u = z{u) 






where 
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minimizes the statistical variance of the free energy estimation for u = AGab- 

The right-hand side of equation ([3]) is expected to provide a constant function near 
the solution of the transcendental equation, as shown in Figs. 2c and 3c for each branch. 

Free energy recovery at zero force 

The solution of the Bennett acceptance ratio method gives the free energy difference 
between state A at Ao and state B at Ai- In order to recover the free energy at zero 
force of each structure with respect to the random coil state, AGo, we need to subtract 
the free energy of stretching the ssDNA, W ss dnA) the free energy of orientation of the 
hairpins stem, VF s tem, and the reversible work performed to stretch the handles and 
displace the bead in the optical trap, Whb' 

AG°ab = &Gab - W ssD na - W stem - W hh (5) 

To compute the work needed to reversibly stretch the ssDNA, W ss dna = / fdx, we 
use the non-extensible worm-like chain elastic model with a persistence length equal to 
1.3±0.2 nm and a contour length equal to the contour length of the molecule. The free 
energy of the stem orientation is evaluated using the freely-jointed chain with a Kuhn 
length equal to the diameter of the hairpin at zero force (Supplementary Section SI) [2]. 
The short length of the handles allows us to estimate the free energy of the handles and 
the bead by integrating a linear FDC along the folded branch from the minimum force 
at Ao to the maximum force at Ai, that is, W^b = (/max — fmm)/2k e E, where k e g is the 
slope of the FDC 0. 



Free energy prediction 

In order to obtain the most stable structure of the DNA molecules under study we use 
the mfold web server [3] . To predict the free energy of formation of each structure we use 
the nearest-neighbor model (Supplementary Section SI). The base pairing free energies 
have been derived in thermal denaturation experiments |3] and independently verified 
in single molecule experiments |18j . 
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Fig. Legends 

Fig. 1. 
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Schematic illustration of the conflgurational space, molecular free energy 
landscapes and experimental setup, a) Conflgurational space partitioned into re- 
gions that correspond to different molecular kinetic states. Inside a partially equilibrated 
region (N, I, U, or M) configurations are sampled according to the Boltzmamr-Gibbs 
distribution; in contrast, the statistical weights of the regions do not necessarily follow 
an equilibrium distribution, b) Schematic free energy landscapes for DNA sequences ex- 
hibiting an intermediate kinetic state on-pathway (II, top) and a misfolded kinetic state 
off-pathway (Ml, bottom). For Ml, the free energy landscapes of the native (black) and 
the misfolded (gray) structures are sketched. Free energies and extensions are shown 
in arbitrary units, c) Experimental setup (not to scale). One bead is immobilized in a 
micropipette by air suction, while the other is captured in an optical trap produced by 
a highly stable dual-beam optical tweezers apparatus [18]. 
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Hairpin II, with an intermediate state, a) Schematic illustration of the mechan- 
ical unfolding/folding pathway (top) and collection of unfolding (red) and folding (blue) 
FDCs recorded at 60 nm/s exhibiting three branches of force corresponding to three 
states: N, I and U. Twelve molecules were pulled at 60 nm/s; between 80 and 385 cycles 
were collected for each molecule. Inset: The gray area indicates the work measured for a 
given FDC between Ao=0 and Ai=55.6 nm. b) Partial work histograms for work values 
measured between Ao=0, where j4=N, and Ai=55.6 nm, where i?=N (red, top panel), 
B=I (green, middle panel) or B=XJ (blue, bottom panel). Dark colors refer to unfolding 
work distributions and light colors to folding work distributions, c) The acceptance ratio 
method (Methods) applied to the work measurements shown in panel b to obtain the free 
energy differences AGnb for B=N (red squares), I (green circles) and U (blue triangles). 
d) Reconstruction of the free energy branches for states N, I and U (color code as in c) 
obtained by fixing Aq and letting Ai change between 45 and 65 nm. The black-dashed 



curve is the full free energy of the system, AG 



inset: Contribution of the prefactor log(^ _s> / < Pr B ) t° the free energy of each kinetic 
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state as a function of A. Right inset: Difference between the free energy of each state, 
AGnb, and the full equilibrium free energy, AG. Error bars in panels b, c, d indicate 
the standard statistical deviation computed over 385 cycles for a given molecule. These 
were obtained using the bootstrap method. 

Fig. 3. 
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Hairpin Ml, with one misfolded state, a) Schematic illustration of the me- 
chanical unfolding/folding pathways (top) and patterns identified in the FDC. For the 
sake of clarity curves are shifted. Left: unfolding (purple) and folding (green) FDC 
for the N state show a force rip (~ 0.5 pN) around 15 pN. Right: unfolding (red) and 
folding (blue) curves are the FDC for the M state. In the analysis we used data from 
5 molecules, pulled at 40, 90 and 125 nm/s. The number of unfolding/folding cycles 
performed varies between 80 and 270 for each dataset. Inset: Unfolding curves for N 
and M without shift (data filtered). At low forces M has larger extension than N. b) 
Partial work histograms for work values measured at 125 nm/s between Ao=0, where 
A=N (red, top panel) or M (blue, bottom panel) and Ai=130 nm, where B=XJ. Dark 
colors refer to unfolding and light colors to folding work distributions, c) The accep- 
tance ratio method is applied to work values measured at different pulling speeds (40, 90 
and 125 nm/s) to extract AGnu an d AGmu- Different colors refer to results obtained 
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for different molecules pulled at different speeds, d) Histograms of the free energy of 
formation of N (red) and M (blue) obtained for different molecules. Error bars in panels 
b, d indicate the standard statistical deviation computed over 270 cycles for a given 
molecule. These were obtained using the bootstrap method. 
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Hairpin 12, with two intermediate states, a) Schematic illustration of the 
mechanical unfolding/folding pathway (top) and collection of unfolding (red) and folding 
(blue) FDCs recorded at 60 nm/s: the four branches of force correspond to the four 
states N, F, I" and U. Curves plotted using light colors are measured by setting extreme 
values of A such that the molecule equilibrates at N and U at the initial and final 
pulling conditions respectively. In curves plotted using dark colors Ai was chosen so 
that at the end of the F process and at the beginning of the R process the molecule 
can be found in any state: N (two patterns are observed in the reverse trajectories with 
different degrees of dissipation), I', I" or U. Curves corresponding to different transitions 
are shifted for the sake of clarity. 7 molecules were pulled at 60 nm/s; between 100 
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and 400 cycles were collected for each molecule, b) Partial work histograms for work 
values measured between Ao=0, where j4=N, and Ai=183 nm, where S=N (red, first 
panel), B=Y (green, second panel), B=V (blue, third panel) or B=XJ (purple, bottom 
panel). Dark colors refer to forward work distributions and light colors to reversed work 
distributions, c) Free energy branches of states N (red squares), I' (green circles), I" 
(blue triangles) and U (purple diamonds) measured relative to the full free energy of the 



system, AG = -k B T log^2 B= 



:N,I',I",U ' 



NB ' B . Error bars in panels b, d indicate 



the standard statistical deviation computed over 400 cycles for a given molecule. These 
were obtained using the bootstrap method. 
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Hairpin M2, with two misfolded states, a) Patterns identified in the FDC 
with corresponding unfolding/ folding pathways. N and M" share a piece of hairpin in 
their folded conformation (orange). Left: unfolding (red) and folding (blue) FDCs for 
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state N show a force rip (~2 pN) around 15 pN. Middle: unfolding (purple) and folding 
(green) FDCs for state M' show no hysteresis. Right: folding (gray) FDCs for state 
M" are identical to the ones measured for M', whereas unfolding (cyan) FDCs show a 
rescue to the N state. Six molecules were pulled, obtaining a minimum of 40 cycles and 
a maximum of 100 at 60 nm/s. b) Partial work histograms for work values measured 
between Ao=0, where A=N (red, top panel), M' (blue, middle panel) or M" (green, 
bottom panel), and Ai=230 nm, where B=XJ. Dark colors refer to unfolding and light 
colors to folding work distributions, c) Free energy branches of states N (red squares), 
M' (blue circles), M" (green triangles) and U (purple diamonds) measured relative to 
the full free energy of the system, AG = — kBTlogY^A=T<i M' M" u e~ AGAU ' kBT . Error 
bars in panels b, c indicate the standard statistical computed over 100 cycles for a given 
molecule. These were obtained using the bootstrap method. 
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Tables 







EFR 


Hopping 


UO 


Unzipping 


11 


AG NU 


30±3 
55±3 


31±2 
61±2 


30.5 
60.1 


27.1 
56.1 


Ml 


AG NU 


47±2 
62±3 


46±3 

58±3 


49.6 
60.2 


46.9 
57.2 


12 


AG NI , 

AG NI » 
AG NU 


40±6 
80±7 

125±7 


- 


41.9 
83.8 
138.0 


39.4 

78.8 
129.1 


M2 


AG NU 

AG M //u 


94±2 
60±3 
70±3 


- 


92.9 
62.0 
72.3 


87.4 
57.4 
67.9 



Table 1: Free energies of formation evaluated using the EFR and equilibrium-based 
methods (Supplementary Section S3) for the DNA sequences studied in this paper com- 
pared to predictions based on the nearest-neighbor (NN) model using the UO set of 
parameters 0] and data obtained from unzipping experiments |18j . Error bars in 
the first and second columns contain the standard deviations over different molecules 
(statistics given in figure captions) and systematic errors in the calibration of the in- 
strument (5% in force and distance). Among the two, the greatest contribution turns 
out to be the force and distance calibration errors. These errors stem from appropriate 
conversion factors between measured voltages in the detectors of forces and distances 
and are multiplicative (Section S2 in [18J). This means that our free energy numbers 
have an absolute value as currently indicated here, and the error in the ratio between ex- 
perimental free energies is smaller. Such error bars are compatible with the discrepancy 
observed between the different predictions provided by the NN model. An extension of 
this table is given in the Supplementary Section S8. 
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Supplementary Information 

SI Evaluation of Free energy landscapes 

Each molecular structure (native or misfolded) has a fixed number N of basepairs in 
the folded state. Along the unfolding pathway, the molecule can explore multiple con- 
formations; their number grows exponentially with the number of bases, and so grows 
the number of potentially stable intermediate states. To model free energy landscapes 
(FEL) only sequential configurations are taken into consideration [lj. Each configuration 
is labeled by the number of open basepairs n: n=0 corresponds to the folded state (all 
basepairs are formed) and n = N to the fully unzipped unfolded state. 
For a given value n and a force / the free energy is given by [SJ: 

AG n (/) = AG° + AG- DNA (/) + AG d n (f) (S6a) 

AG ssDNA (/) = / F^ NA (x')dx> - f Xn (f) 

Jo 

x n (f')df (S6b) 





r x dU) 
AG d n (f) = / F d (x')dx' - fx d (f) 

Jo 



x d {f)df (S6c) 

o 

where AG^ is the free energy of formation of the n th configuration at zero force that 
can be estimated using the nearest-neighbor (NN) model and the unified-oligonucleotide 
(UO) set of parameters derived from bulk experiments [31 3] or data derived from single- 
molecule unzipping experiments [5J. AG^f (/) is the elastic free energy at force / 
of the released ssDNA for the configuration n, x n (f) being its equilibrium end-to-end 
distance projected along the force axis and -F^dna^) ^ ne e Q ua tion of state of a ssDNA 
polymer of contour length l n . We use the inextensible worm- like chain (WLC) model 
with a persistence length P = 1.3 ± 0.2 nm and an inter-phosphate distance equal to 
0.59 nm/base [51 El El [9] . AG^(/) is the energetic contribution due to the orientation 
of the hairpin double helix (modeled as a single dipole of length d) along the force axis. 
Here we use an inextensible freely-jointed chain (FJC) model, with Kuhn and contour 
lengths equal to the diameter of the double helix, d ~2.0 nm 0[8]. 

If the structure under consideration is made of more than one hairpin (state M for 
Ml; N, I' and I" for 12; or M' and M" for M2) multiple sequential configurations can 
take place for an intermediate value of n. In this case, an exponential Boltzmann sum 
of the different contributions gives the free energy AG n (/). 

Following, the FEL of the different folded structures under consideration are shown. 
II and 12 fold only into one native structure and intermediates are revealed along the FEL 
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(Figs. SI and S2). Molecule Ml can fold via two different pathways into two unrelated 



structures: the native (N, Fig. |S3[) or the misfolded (M, Fig. S4). Molecule M2 can fold 



Fig. S6 or M" in Fig. S7). 



into one native structure (N in Fig. S5), and two different misfolded structures (M' in 



S2 Free energy recovery from equilibrium experiments 

S2.1 Molecule II 

In the equilibrium-based hopping experiments the trap-pipette distance A is kept sta- 
tionary and the molecule executes transitions between the three different states (N, I, 



U). A typical trace is shown in Fig. S8 



i.. 



The histogram of the measured force along an equilibrium trace can be fitted to a 
sum of three Gaussians: 



P(f) 



w N 



2na 2 N 



i (/-(/jy)r 
-e n 



+ 



Wj 



2ira 2 j 



1 (f-ifi)r 

~2 ^ 

-e i 



+ 



wu 



2-nal 



i u-(fu)r 
-e v 



(ST) 



where u>n, wi and wjj are the relative weights of each state N, I and U respectively 
(wn + wj + wu = 1); (/)jv, (/)/ and (f)u are their average forces; and cr/v, 07 and o\j 
are the standard deviations. The free energy of formation of I and U (relative to N) 
can be obtained using the detailed balance condition, by relating the logarithm of the 
relative weights of the Gaussians to the free energy differences between states. To get 
the free energy at zero force we need to subtract the elastic contribution of each state 
(derived in section S2) at the average force of the equilibrium trace (/}: 



AG 



NI 



-k B T log 



W N 



■(f) 

~ I xj(f)df 

'0 



AG NU 



wnJ 



x d (f)df+ / Xu (f)df 
o Jo 



(S8) 



(S9) 



Results, summarized in Table 1 in the main paper, are in agreement with free energy 
predictions using the NN model with the UO set of parameters and data from unzipping 

lasusj. 



S2.2 Molecule Ml 

In the case of equilibrium experiments with Ml, kinetics are very fast and some tran- 
sitions are missed. Hopping traces in Ml reveal two clearly separated hopping regions 



(Fig. S9): in one region states N and an intermediate In coexist; in the other region 
states M, U and another intermediate Im coexist. A closer look of the traces at the 
interphase between both regions shows that there are no trajectories directly connecting 



NandM (Fig. S9). 
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The histogram of the measured force along an equilibrium trace can be fitted to a 
sum of five Gaussians: 

2 j (J-(fl N )) (f-thr)) 2 



l (/-(/iV>r ,., -A " "JN" l (/-(/j7» 
"2 ?2 W/^ 2 „3 WJJ ~2 ?2 

1 (f - {f lM )] 1 (/-</jv>) 2 

H e 7 m -| e iv (S10) 

V 2vrcJ L v 2vrcJ ^ 

where to*, i =N, Ijv, U, Im, M are the relative weights of each state (^ Wi = l); (f)i are 
their average forces; and ctj are the standard deviations. The free energy of formation of N 
and M (relative to U) can be obtained using the detailed balance condition, by relating 
the logarithm of the relative weights of the Gaussians to the free energy differences 
between states. To get the free energy at zero force we need to subtract the elastic 
contribution of each state at the average force of the equilibrium trace (/) : 

AG° NV = -k B T log ( ^-) - [ x d (f)df+ [ xu(f)df (Sll) 

V^n/ Jo Jo 

( V1tt\ fW fW fW 

AG° MlJ = -k B Tlog -^ - 2 / x d (f)df - / x M (f)df + / Xu (f)df (S12) 
\wmJ Jo Jo Jo 

Results, summarized in Table 1 in the main paper, are in agreement with free energy 
predictions using the NN model with the UO set of parameters and data from unzipping 
[31 IU [5]. Intermediate states are not characterized. 

S2.3 Molecule 12 

Equilibrium-based hopping experiments reveal four levels of force, corresponding to the 
four different states N, I', I" and U. Hopping between all four conformations was never 
observed in the accessible experimental time scales: at high forces hopping between I', 
I" and U (but not N) occurs (Fig. SIO^,); at low forces once the molecule reaches N it 



never leaves that conformation precluding equilibrium based free energy measurements 
(FiglSlOfc). 



In the absence of coexistence between the four states the free energy of formation of I', 
I" and U (relative to N) was not evaluated using equilibrium-based hopping experiments. 

S2.4 Molecule M2 

The network of intermediate states present in the folding pathways of N, M' and M" 



taking into account only sequential configurations is extremely complex (Fig. Sll). In 



fact, the molecule can change state from M' or M" to N without going through U. 
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Figure SI: Molecule II. Sequence, free energy landscape and FDC. a. Sequence 
and secondary structure. The native state of II is characterized by the presence of an 
internal loop in the middle of the stem that favors the existence of an intermediate (I) 
along the unfolding pathway of the hairpin, b. Free energy landscape evaluated at the 
coexistence force (where the folded state, n = 0, and the unfolded state, n = N, have 
equal free energies). I contains four different configurations, n = 11 — 14, with the same 
free energy AG n (f c ). c. Example of unfolding FDC (red) and refolding FDC (blue, 
yellow). The branches of force are indicated with black straight state-lines. Each branch 
is assigned to one state of the molecule: N, I or U. Inset: Detail of an unfolding curve 
where some points (indicated with red arrows) can be assigned either to I or U depending 
on the size of the running average used for the classification of the data points. 
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Figure S2: Molecule 12. Sequence, free energy landscape and FDC. a. Sequence 
and secondary structure. The native state of 12 is characterized by the presence of a three 
way junction that favors the existence of two intermediates (I' and I") along the unfolding 
pathway of the hairpin, b. Free energy landscape calculated at the coexistence force 
(15.08 pN) predicted by the NN model using the UO parameters. Two intermediates, 
hereafter referred as I' and I" , appear along the unfolding pathway surrounded by high 
kinetic barriers (10 k-oT). V is found at n = 21 and corresponds to a structure where 
all the stem is open and the two hairpins of the bifurcation are closed. I" is found 
at n = 37 and corresponds to two identical structures, where all the stem and one 
hairpin in the bifurcation are open. As this two configurations have identical molecular 
extension, we will consider I" as a single state. Blue squares are the free energies of all 
the possible sequential configurations that have a given number of n of unzipped base 
pairs distributed between the two hairpins. Red circles are the Boltzmann average (i.e. 
the mean free energy potential) taken by summing over all configurations constrained by 
a given number of unzipped base pairs n. c. Example of unfolding (red) and refolding 
(blue) FDC along one trajectory. The four branches of force are assigned to the four 
states N, I', I" and U. 
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Figure S3: Molecule Ml, N state. Sequence, free energy landscape and FDC. 

a. Sequence and secondary structure of N, which consists of a single-stem hairpin, b. 
Free energy landscape evaluated at the coexistence force (14.75 pN, red squares), at 
16 pN and at 14 pN. At the coexistence force an intermediate located at n = 13 with 
a very low kinetic barrier to the native state (n = 0) can be seen. This intermediate 
may play a role in both unfolding and folding FDC (panel c), where force fluctuations 
increase in the vicinity of unfolding and folding rupture forces. For the sake of simplicity 
this intermediate is not characterized in our study, c. Example of unfolding (red) and 
refolding (blue) FDC. 
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Figure S4: Molecule Ml, M state. Sequence, free energy landscape and FDC. 

a. Sequence and secondary structure of M, which consists of two hairpins serially con- 
nected by seven thymines, b. Free energy landscape evaluated at the coexistence force. 
Blue squares are the free energies of all the possible sequential configurations of the 
structure. Red circles are the exponential Boltzmann average over configurations con- 
strained by a given number of unzipped basepairs n. An intermediate state surrounded 
by high kinetic barriers (~8 k^T) appears along the unfolding pathway at n = 10. Con- 
figurations that dominate the exponential sum in the Boltzmann average for this value 
of n are (10,0) and (0,10), where one of the hairpins is unzipped and the other is folded. 
For the sake of simplicity this intermediate is not characterized, c. Example of unfolding 
(red) and refolding (blue) FDC. 
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Figure S5: Molecule M2, N state. Sequence, free energy landscape and FDC. 

a. Sequence and secondary structure of N. b. Free energy landscape evaluated at the 
coexistence force (13.52 pN), 16 pN and 11 pN. c. Example of unfolding (red) and 
refolding (blue) FDC. 
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Figure S6: Molecule M2, M' state. Sequence, free energy landscape and FDC. 

a. Sequence and secondary structure of M', which consists of four small hairpins that do 
not share any structure with N . b. Free energy landscape evaluated at the coexistence 
force. Blue squares are the free energies of all the possible sequential configurations 
of the structure. Red circles are the Boltzmann average taken over all configurations 
constrained by a given number of unzipped base pairs n. Several intermediates separated 
by low kinetic barriers (~5 k-gT) appear: one is located at n = 6, being the most stable 
configuration (0,6,0,0) where the first, third and fourth hairpins are folded and the second 
hairpin is fully unfolded. Other intermediates are found at n = 12 (where (0,6,6,0) and 
(0,0,12,0) dominate the Boltzmann average), n = 18 (0,6,12,0) and n = 25 ((7,6,12,0) 
and (0,6,12,7)). c. Example of unfolding (red) and refolding (blue) FDC. The presence 
of many transition states with low kinetic barriers implies an almost reversible pattern 
for the FDC. 
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Figure S7: Molecule M2, M" state. Sequence, free energy landscape and FDC. 

a. Sequence and secondary structure of M" , which consists of three hairpins, one hav- 
ing 21 basepairs in common with N (Fig. S5). b. Free energy landscape evaluated at 
the coexistence force. Blue squares are the free energies of all the possible sequential 
configurations of the structure. Red circles are the Boltzmann average taken over all 
configurations constrained by a given number of unzipped base pairs n. The FEL evalu- 
ated at the coexistence force reveals two intermediates at n = 22 (configuration (0,22,0)) 
and n = 29 (configurations (7,22,0) and (0,22,7)). c. Example of unfolding (red) and 
refolding (blue) FDC. 



26 



£ 

a 




10 
Time (s) 



20 25 0.4 Q.h 

P(f) 



Figure S8: Hopping experiments with hairpin II. In these experiments the trap- 
pipette distance A is kept stationary, a. Sample trace during 25 seconds. In red we 
show data at full resolution and in black we show data filtered. The molecule executes 
transitions between N (at ~14.5 pN), I (at ~13.8 pN) and U (at ~13 pN). b. Histogram 
of the measured forces during the hopping experiment. The free energy of formation of 
each state can be obtained by fitting the histogram to a sum of three Gaussians and 
using the detailed balance condition, by relating the logarithm of the relative weights of 
the Gaussians to the free energy differences between the states. 
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Figure S9: Hopping experiments 'with hairpin Ml. In these experiments the trap- 
pipette distance A is kept stationary, a. Sample trace at full resolution. The molecule 
executes transitions between N (~14.1 pN), M (~13.9 pN), two intermediate states In 
and Im on-pathway to the native and to the misfolded state respectively (~13.5 pN in 
both cases since the molecular extension is identical) and U (~13.0 pN). c. Detail of 
the boxed region in panel a. A closer look of the traces at the interphase between both 
regions shows that there are not trajectories directly connecting N and M. b. Histogram 
of the measured forces during the hopping experiment. The free energy of formation of 
each state can be obtained using a fit to a sum of five Gaussians. 
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Figure S10: Pulling and hopping experiments with hairpin 12. a. Example 
of a passive hopping trace measured at high forces: the molecule is initially set to N. 
Once the molecule partially unfolds (after ~15 in the figure), the system subsequently 
hops between states I', I" and U, but never folds back to N. b. Example of a hopping 
experiment at low forces: The molecule is initially set to I' (it was not possible to stabilize 
the molecule in neither U nor I" states at such low forces) and once the molecule folds 
to N it never escapes out. Equilibrium hopping experiments were repeated at different 
forces and the native state was never observed to coexist with any other state. 
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Figure Sll: Folding pathways for N, M' and M". a) Folded structures for M2 
according to Mfold [3, 4J: M' is made by four small hairpins hereafter referred as Hi, 
H2, H3 and H4; M" is made by hairpins H3 and H4 (also present in M'), and hairpin 
H5; and N is made by a single and long hairpin that contains H5. b) The folding 
pathway to each state can be modeled by the sequential folding of Hi, H2, H3, H4 or 
H5. Starting from U a first hairpin may form (Hj, i=l,2,3,4 or 5). Depending on the 
first hairpin, a second hairpin may form (Hjj=Hj+Hj). Folding proceeds to the next 
level of structures (Rijk) down to the final states M', M", N. Black (red) arrows denote 
pathways that need to overcome a kinetic barrier of ~5 (~10) k^T. The blue arrow 
denotes the transition N^Hs , which might be mediated by an intermediate on-pathway 
predicted in the free energy landscape, c) Free energy landscape along the different 
folding pathways calculated at a constant force of 11 pN using the NN model and the 
UO set of parameters. Folding pathways with barriers that are higher than 10 k-gT are 
plotted in red, and the transition towards N is plotted in blue. Green circles denote 
first-level configurations of the type Hj, blue triangles denote second- level configurations 
of the type H^ , purple diamonds denote third-level configurations of the type H^ and 
black squares indicate the initial state U and the final states N, M' and M". 
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In the case of M2 kinetics in equilibrium experiments are not very fast (Fig. S12) 



but it is very difficult to identify different conformations and unfolding/folding pathways 
along a hopping trace (due to the nearly identical molecular extensions of some states) . 



In addition, as shown in Fig. S12, the complexity of the network of states for M2 is such 
that there are pathways connecting N, M' and M" that do not go through U. 

Therefore, it is not possible to measure the free energies of N, M' and M" (relative 
to U) and to unravel unfolding pathways from equilibrium-based hopping experiments. 
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Figure S12: Hopping experiments for hairpin M2. At an intermediate range of 
forces (10-13pN) the hairpin shows hopping between N (mean force 12.7 pN), an inter- 
mediate state (mean force 11.7 pN) and a mixture several states ( mean force llpN). 
The latter mixture is potentially composed of U, M', M" and other intermediates (boxed 
region, see also Fig. |S11[ ) that exhibit nearly the same molecular extension being very 
difficult to distinguish (rectangular box). 



S3 Identification of folding/unfolding pathways along FDCs 



In order to prove that molecules fold into the structures summarized in Figs. |Sl||S7| we 
measure the molecular extension released/absorbed, Ax m , along FDC as: 



AXr 



A/ 

&eff' 



(S13) 



where A/ is the change in force along the transition and k e g is the effective stiffness 
of the molecular construct (dsDNA handles plus optical trap), measured as the slope 
of the FDC before the transition occurs. In Tables [S2]|S5| the experimentally measured 



Ax m for the different molecules are compared to predictions obtained using the worm- 
like-chain (WLC) model with a persistence length equal to 1.35 nm. Results match 
theoretical predictions. 
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A/ ( P N) Ax^ p (nm) Ax% LV (nm) 
N->1 0.76±0.04 12±1 rXo 

1^1 0.57±0.03 11±1 11.2 

Table S2: Change in molecular extension measured for different transitions in 
hairpin II. Results are in agreement with predictions obtained using the WLC model. 

Statistical errors are insignificant compared to systematic error, which we take equal to 

5%. 



A/ (pN) Axg p (nm) Ax% LV (nnij~ 
N^U 1.20±0.06 24±1 23.0 

M->U (*) 0.80±0.08 16±2 18.1 

Table S3: Change in molecular extension measured for different transitions 
in hairpin Ml. Results are in agreement with predictions obtained using the WLC 
model. Statistical errors are insignificant compared to systematic error, which we take 
equal to 5%. (*) FDCs do not show a sudden jump in force for the transition M— >XJ. 
(A/) was then calculated from the relative shift between each force branch M and U 
(see the inset of Fig. 3a at low forces in main text). 



A/ (pN) Ax^ p (nm) Ax% LLj (nm) 

N->T 1.30±0.07 20±1 2L1 

P-»F 0.82±0.05 13±1 14.7 

I"->U 0.80±0.05 14±1 14.7 

Table S4: Change in molecular extension measured for different transitions in 
hairpin 12. Results are in agreement with predictions obtained using the WLC model. 
Statistical errors are insignificant compared to systematic error, which we take equal to 

5%. 





A/ (pN) 


Axm P (nm) 


Ax« LU (nm) 


N^U 


1.50±0.08 


35±2 


38.5 


M"^N 


0.80±0.05 


18±1 


15.9 


M^U (*) 


1.3±0.2 


28±4 


30.8 



Table S5: Change in molecular extension measured for different transitions 
in hairpin M2. Results are in agreement with predictions obtained using the WLC 
model. Statistical errors are insignificant compared to systematic error, which we take 
equal to 5%. (*) FDCs do not show a sudden jump in force for the transition M— >XJ. 
(A/) was calculated from the relative shift between each force branch M and U (as in 
Ml), and values of for M' and for M". 
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S4 Validation of EFR 



From the extended fluctuation relation (EFR) [10] we can write: 

\W) \ \Y AGab 



log 



pA^B{ 
r F 



P^ B (-W)J k B T k B T 



log 






(S14) 



If we plot log 



pA-+B 






(W) 



<(-W) 



as a function of the work in k B T units we obtain a straight 



line with slope equal to 1 as a direct proof of the validity of the EFR. Fig. S13 show 
experimental tests of the validity of the EFR for molecule II taking ^4=N and B=N, I or 
U. The partial work distributions used are shown in Fig. 2b in the main paper. Within 
experimental errors the EFR holds for the three different kinetic states. 

a b c 




154 155 156 
W (k B T) 



50 151 152 153 154 155 156 157 

W {k B T) 




153 154 155 156 

W (k B T) 



157 



Figure S13: EFR test for II. Validation of the EFR using the partial work distributions 
shown in Fig. 2 for states N (a), I (b) and U (c). Error bars were obtained using the 
bootstrap method. 



S5 Potential sources of error: classification of states along 
FDCs 

S5.1 Molecules with intermediates on-pathway 



Measured FDC for molecules II and 12 show different force branches (Figs. SI z and S2 



Each branch corresponds to one state. By representing each branch with a characteristic 
straight state-line (black lines in Fig. SI z) it is possible to assign a state to each measured 



data point (A, /). This is done by determining the nearest state-line to each measured 
point. The drawback of this method stems from large force fluctuations (Fig. SI ;, inset): 



due to the finite acquisition rate, 1 kHz, sometimes there are isolated points that can be 
assigned to the wrong state. As a consequence, an error can be introduced in the fraction 
4>p~^ B and in the partial work distributions Pp~^ B (W) and P^ B {— W), where ^4=N 
and JB=N, I or U (</>/j = 1 since at the chosen value of Ao all reversed trajectories end 
in state N). 
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In order to study this effect, different sizes of a running average (1, 5 or 10 points) of 



experimental data for molecule II are considered. Results for the prefactor log (<ftp /4>r'~' ' 



evaluated taking Aq=0 nm (A=N) and varying Ai between 20 and 80 nm can be seen in 



Fig. S14i for B=H (red), I (black) and U (blue). A convergence is observed for running 
averages between five and ten points. The error made in the evaluation of free energies 
can be large for intermediate states under short lifetime conditions. For instance, in the 
case of II an error of the order of 2 k-gT is made only for I at low values of \\. 
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Figure S14: Molecule II. 



cation, a. The quantity log</>^ _ShS 



Effect of running boxcar averages on state classifi- 

^~ B as a function of Ai is evaluated for different 
sizes of the running average of experimental raw data: 1 point (asterisks), 5 points (open 
squares) and 10 points (solid circles). B = N for red, B = I for black and B = U for 
blue symbols, b. Difference between the free energy of each state AGnb and the full 
free energy, AG = —k^T log J2 B= ^ lv e~^ GNB ' kBT , evaluated at different values of Ai 
for different sizes of the running average. Color code as in a. 



S5.2 Molecules with misfolded states off-pathway 

From the whole FDC pattern we can identify if the hairpin has folded into states N 
or M. Therefore, no error is made when evaluating 4>r~ (4'f~^ B= ^ always) for a high 
value of Ai (where molecule is unfolded) and a low value of Ao (where molecule is folded, 
either in N or M). 

In the case of M2 we assume that unfolding curves starting from a misfolded state 
and showing rescue to N around 9.5 pN (Fig. S7;) start in M" [TT]. Folding curves 



preceding unfolding curves that are rescued by N are assigned to fold into M"; and 
folding curves preceding unfolding curves without rescue are assigned to fold into M'. 
A crucial question to justify this assumption is how far from equilibrium the system is: 
if we are under quasi-static conditions the molecule has more time to explore the free 
energy landscape and to overcome high kinetic barriers, preferentially folding into N. 
On the other hand, if we are far from equilibrium the molecule will mostly fold into M' 
because its folding pathway encounters lower kinetic barriers. In fact, as shown in the 



folding network of M2 sketched in Fig. Sll the molecule can change conformation from 
M" to N via two intermediates (H45 or H35 and H5) without crossing too high kinetic 
barriers (~5 k^T, shown as black arrows), whereas the transition from M' to N involves 
several intermediates (H-ijk, Hjj and Hj) and at least one high kinetic barrier (~10 k^T, 
shown as red arrows). Consequently, at an intermediate value of the pulling speed (not 
too far from equilibrium), unfolding events from a misfolded state that show rescue to 
N probably correspond to the unfolding of M". To conclude, it is crucial to be out of 
equilibrium to favor misfolded states, but not too far so that rescue to N from M" is 
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preferential. 

If the distinction between M' and M" is not made and a coarse grained state M 
is considered (M=M'UM"), a free energy for M that is nearly equal to that of M' is 
obtained. That is in accordance with the fact that the Boltzmann average of the free 
energies of M' and M" , 



AG 



M 



-k B T log 



exp 



AGm' 
k B T 



+ exp 



AG 



M>' 



k B T 



AG M ' = 61.97fc B T 



(S15) 



is dominated by the lowest free energy among the two. In fact, AGm" exceeds AGm' 
by 10 k B T units (Table 1, main text). 



S6 Effect of neglecting the term log ((f) p /(f) 



-B" 



R 



The consequences of neglecting the correction log (</> I p - ^/</ >I fj' <-S ) can be seen in Fig. 



S15 1, where results for the acceptance ratio method applied to work measurements 
between Ao = and Ai = 55.6 nm are shown. Note that the relative stability of the 
three states N, I and U change as compared to Fig. 2c in the main paper. The same trend 
can be observed in Fig. S15b , where the free energy branches of the three kinetic states 
calculated by neglecting the correction are shown. These results should be compared 
with Fig. 2d (right inset) in the main paper. 




146 150 154 158 162 
u (k B T) 



E-H 

< 2 







20 



! N / 


- 


i i i 





40 60 

Ai (nm) 



NO 



Figure S15: Molecule II. Effect of neglecting the prefactor 0p^ B /0R^ B . a. 

The acceptance ratio method is applied to the work values measured between Ao = to 
Ai = 55.6 nm. b. Difference between the free energy of each state, AGnb, B = N (red), 
I (black), U (blue) and the free energy of the system AG. Note the change with respect 
to the results shown in Fig. 2c and 2d (right panel) in the main text. Error bars were 
obtained from the bootstrap method. 
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S7 Contributions to the evaluation of the free energy at 
zero force 







11 






Ml 




A =0 


, Ai=55.6 nm 


A =0, Ai=130.0 nm 


N->N 


N-)-I 


N^U 


N^U 


M-^U 


(W)£+ B 


156±1 


155±1 


154±1 


468±2 


260±1 


(W)^ B 


155±1 


152±1 


149±1 


465±1 


459±1 


log {<t>p B /<j>jr B ) 


-0.11 


-2.77 


-3.12 


2.16 


0.12 


AG™ 


156±1 


157±1 


156±1 


464±2 


459±1 


AG Cmoks 


156±1 


154±1 


152±1 


466±2 


460±1 


(w diss )^ B 


Oil 


2±1 


1±1 


3±2 


Oil 


(w diss )i^ B 


Oil 


5±1 


6±1 


1±2 


Oil 


w hb 


156±3 


115±3 


83±3 


381±3 


392±3 


WssBNA 





9.3 


18.4 


22.4 


21.0 


Ar ,0(EFR) 


0±3 


30±3 


55±3 


61±3 


49±3 
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M2 






A 


o=0, Ai= 


= 148.0 nm 


A =0 


, Ai=230.0 nm 


N-^N 


N-*r 


N-»I" 


N^U 


N^U 


M'^U 


M"^U 


{W)jr+» 


544±1 


543±1 


541±1 


544±1 


623±1 


583±1 


611±1 


{W)f B 


533±3 


508±1 


507±1 


503±1 


594±1 


579±1 


580±1 


log {^ B /^ B ) 


-0.18 


-3.03 


-2.42 


-3.72 


0.71 


1.75 


1.1 


AG™ 


543±3 


527±1 


526±1 


528±1 


612±1 


579±1 


590±1 


AG Crooks 


544±3 


525±1 


525±1 


526±1 


613±1 


581±1 


592±1 


(W diss )^ B 


1±1 


6±1 


15±1 


16±1 


11±1 


4±1 


21±1 


(w diss )£*- B 


10±3 


19±1 


21±1 


25±1 


18±1 


Oil 


10±1 


w hb 


543±4 


469±6 


417±6 


365±6 


480±2 


487±3 


485±3 


W^ssDNA 





18.9 


30.5 


41.7 


37.5 


31.5 


35.5 


ArY 0(EFR) 


0±4 


40±6 


80±7 


125±7 


95±3 


60±3 


69±3 



Table S6: Relevant energy contributions. For each molecule, we show the measured 
average forward and reverse works {W)p B and (W)^~ B , the contribution to the free 
energy introduced by the correction term log {4>f* B /4>ii~ B )i the free energy of the 
system obtained using the EFR AG™ (equation 1) and the Crooks fluctuation theorem 
(no correction term included) AG9g° [H], the average forward and reversed dissipated 
works (W d i ss )p B and (Wdiss)ij B , the reversible work performed to stretch the handles 
and to displace the bead in the optical trap Whb, an d reversible work needed for the 
stretching of the ssDNA and the orientation of the hairpin stem W ss dna (see Methods). 
Finally, we present the free energy estimation at zero force obtained using the EFR, 



AG' 



0(EFR) 
AB 



Error bars contain statistical and systematic errors. All the magnitudes are 



given in k^T. 
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S8 Extended Jarzynski equality 

If we multiply equation 1 with the reverse partial work distribution Pr*~ (— W) and 
integrate over the work we obtain the Extended Jarzynski equality (EJE): 



ag ab aA-^B I w \A->B 

e ~'- T =5jp (•"'■*), (S16) 

The EJE allows to evaluate AGab using only the forward partial work distribution. 

In the case of 12, the EJE was applied to recover the free energy branches of the 
four states N, I', I" and U. Ao was fixed at a value where A=N for all trajectories, and 
therefore 0jj = 1, and Ai was varied between 100 and 200 nm. The free energy of 
each state at a given value of Ai can be evaluated using: 



AGnb = -k B T log <fir>» - k B T log ( e V \ ( S17 ) 



w x n^b 



Similarly to the Jarzynski equality [13] , the EJE is strongly biased |14| . Using pulling 
experiments where kinetic states are partially equilibrated at Ai (dark curves in Fig. 4a) 
we can evaluate the free energy difference AGnb using the EFR (equation (1)) and 
the EJE (equation (S17|)) for the four different states. The difference between the two 



magnitudes is an estimation of the bias (Table S7). 



B AG*™ (k B T) AG^ B K (k B T) Bias (k B T) 



N 


543±3 


543±3 


0±4 


r 


545±1 


527±1 


18±2 


i" 


539±1 


526±1 


13±2 


u 


547±1 


528±1 


19±2 



Table S7: Estimation of bias for 12. Free energy difference AGnb obtained from the 
EFR and the EJE, and estimation of the bias, for each state i?=N, I', I" and U. 



To evaluate free energy branches we applied equation (S17) to F processes obtained 
from standard pulling experiments (where no intermediates are observed at extreme 
values of A) by fixing Ao=0 nm and varying Ai between 100 and 200 nm. The estimation 
of the bias, taken independent of Ai, is subtracted from the resulting free energies for 
each state. 
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